Action potential as a pressure pulse propagating in the axoplasm 
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We suggest that the propagation of the action potential is driven by a pressure pulse propagating 
in the axoplasm along the axon length. The pressure pulse mechanically activates Na"'' ion channels 
embedded in the axon membrane. This activation initiates the development of a local membrane 
voltage spike, as in the Hodgkin-Huxley model of the action potential. Extracellular Ca^"*" ions 
influxing during the voltage spike trigger a mechanism that amplifies the pressure pulse and therefore 
prevents its viscous decay. The model is able to explain a number of phenomena that are unexplained 
within the Hodgkin-Huxley framework: the Meyer-Overton rule for the effectiveness of anesthetics, 
as well as various mechanical, optical and thermodynamic phenomena accompanying the action 
potential. The model correctly predicts the velocity of propagation of the nerve impulse, including 
its dependence on axon diameter, degree of myelination and temperature. 



I. INTRODUCTION 

The well-known Hodgkin-Huxley theory ofpropaga- 
tion of the nerve impulse in neuronal axons [l| has re- 
cently been the subject of constructive criticism 0, 0, 0, 
[^. As a purely electrical theory, the Hodgkin-Huxley 
theory is unable to explain a number of mechanical 
and thermodynamic phenomena that are observed syn- 
chronously with a propagating action potential. These 
include changes in nerve dimensions and in the normal 
force exerted by the nerve d, SJl, & [lo| re versible 
changes in temperature and heat [111 [l3 . Il3l Q, [lB| and 
changes in fluorescence intensity and anisotropy of lipid 
membrane markers (l6l [iTj . Importantly, within the elec- 
trical framework, it has not been possible to explain the 
famous Meyer-Overton rule [H, [T^] for the effectiveness 
of anesthetics on nerve fibers. The rule states that the 
effectiveness of an anesthetic is linearly related to the 
solubility of that anesthetic in membranes [l^l . The rule 
is valid over 5 orders of magnitude and holds indepen- 
dently of the chemical identity of the anesthetic; it is 
valid for noble gases such as argon and xenon as well as 
for alkanols. No link is found between the reactiveness 
of anesthetics with ion channels, the major actors in the 
Hodgkin-Huxley theory, and the Meyer-Overton rule. 

The Hodgkin-Huxley mechanism of propagation of the 
nerve impulse pj can be divided into two phases. First, 
after an initial depolarization of a membrane segment to 
the excitation level, a local voltage spike develops. Sec- 
ond, ionic currents flowing through the excited membrane 
depolarize adjacent unexcited membrane segments to the 
excitation voltage, causing propagation of the nerve im- 
pulse. The viability of the first phase of the mechanism 
is little doubted. In the voltage clamp method elec- 
trodes are positioned inside and outside of an axon, al- 
lowing the experimenter to change the membrane voltage 
arbitrarily and to measure the resultant membrane volt- 
age and current. Using this method, it is also possible to 
simultaneously depolarize the membrane along the entire 
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length of the axon. If depolarization reaches a certain 
threshold level, one observes a "membrane" action po- 
tential (voltage spike), which develops concurrently along 
the axon length. Hodgkin and Huxley studied this type 
of action potential in detail flj . They showed that the du- 
ration, magnitude and time course of the voltage spike 
are explained by ionic fluxes across the membrane. The 
fluxes arise due to time- and voltage-dependent changes 
in ion-speciflc membrane permeability. There are no ax- 
ial currents during the membrane action potential, which 
simplifies theoretical modeling of the process. 

The second phase of the Hodgkin-Huxley mechanism, 
depolarization of unexcited membrane segments by ionic 
currents straggling in the axial direction along the axon 
length and thus causing impulse propagation, has not 
been shown experimentally. To test the validity of this 
part of the mechanism, one would ideally wish to block 
axial ionic currents at an axonal cross-section and to ob- 
serve whether a propagated action potential is able to 
cross this boundary. To date, such an experiment has 
not been done. Proof of the propagation mechanism has 
been based on numerical calculations. Using their__elec- 
trical framework, Hodgkin and Huxley calculated [1] the 
velocity of the action potential for the squid giant axon. 
Although the calculated velocity agreed well with experi- 
ment, a number of flaws in the formulation of the problem 
were pointed out [5]; e.g., changes in the membrane ca- 
pacitance due to variation in membrane thickness were 
ignored. Hodgkin and Huxley themselves noted [l[ that 
due to non-zero axial currents "the situation is more com- 
plicated in a propagated action potential" compared to a 
membrane action potential, and "it is not practicable to 
solve as it stands" the equation describing the problem. 
To simplify, they assumed that a steadily propagating 
solution exists, in which, additionally, the shape of the 
voltage spike is preserved in time. No analysis of the 
stability of the solution was made. Even with these sim- 
pliflcations, the numerical calculation was in parts "ex- 
cessively tedious" , and " the solution was continued as a 
membrane action potential" . Clearly, the calculation was 
rather complicated and will be expected to have associ- 
ated error. 

In this paper, we present an alternative theory for the 
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propagation of the action potential. We suggest that its 
propagation is driven by a pressure pulse propagating in 
the axoplasm. This model allows a well-understood cal- 
culation of the action potential velocity. For both myeli- 
nated and unmyelinated fibers of various diameters, the 
velocity predicted by this theory agrees very well with 
experiment. The theory is able to quantitatively explain 
the Meyer-Overton rule as well as various mechanical, 
optical and thermodynamic phenomena that accompany 
the action potential. The speculation that density pulses 
may be related to nerve impulses has been discussed by 
various authors since 1912 [SlJSi HI, [H, HI]. Most 
recently, it has been suggested [J, Q that the action po- 
tential is a soliton propagating in the axon membrane. 
While this model appears to account for the abovemen- 
tioned phenomena, as well as for the Meyer-Overton rule, 
it does not seem to explain changes in the propagation 
velocity that are related to axon diameter [26| . 



II. THE HODGKIN-HUXLEY THEORY 
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FIG. 1: Electrical circuit representing membrane. -Rjva ~ 
1/gNa, Rk = l/gK, Ri = l/gi- Rnu and Rk vary with time 
and membrane potential, the other components are constant. 
From P]. 



According to the Hodgkin-Huxley theory [l| , the prop- 
agation of the action potential is driven by electrical cur- 
rents flowing in the axial direction along the length of 
the axon. In the resting state, ion transporters maintain 
transmembrane concentration gradients, pumping in- 
side the cell and Na"*" outside [20| . The resting membrane 
is more permeable to K"*"; therefore, the electrochemical 
equilibrium maintains a negative potential inside the cell. 
Depolarization of the membrane to a threshold level re- 
sults in a rapid and transient opening of voltage-gated 
Na"*" channels, causing the membrane voltage to spike to 
the Nernst equilibrium for Na"*" ions. This is followed by a 
slower activation of voltage-gated K"*" channels, which re- 
polarizes the membrane to the resting potential. The ac- 
tion potential propagates due to axial spread of electrical 
currents from the site of the voltage spike Jj] . These cur- 
rents depolarize adjacent unexcited regions of the mem- 
brane to the excitation voltage. Within this framework, 
the membrane is described as the electrical circuit shown 
in Fig. 1. The current density across the membrane is 
given by: 
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(1) 



where the four terms on the right-hand side of the equa- 
tion give, respectively, the membrane capacity current, 
the current carried by K"*" ions, the current carried by 
Na+ ions, and the "leakage current" (due to Cl^ and 
other ions) per unit area of membrane. Cm is the mem- 
brane capacitance, and gK, 9Na and gi are, respectively, 
the conductances to K+ ions, Na+ ions and the leakage 
currents. Ek, E^a and Ei are, respectively, the equilib- 
rium potentials for K"*" and Na"*" and the leakage currents. 
gNa and gK vary with the membrane potential U and 
time as parameterized from voltage clamp data [l| , while 



Cm, gi, Ek, E^a and Ei are constant. For a "mem- 
brane" action potential, U is the same at each instant 
over the length of the axon. There is no current along 
the axon axis, and / is always zero. The time course of 
the action potential is obtained by solving eq. (1) numer- 
ically, with / = and the initial depolarization condition 
U = Uo- Hodgkin and Huxley showed convincingly [l| 
that the time course of the membrane action potential is 
indeed explained by behavior in g^a and gpc, as param- 
eterized from the voltage clamp data, and the dynamics 
as described by eq. (1). 

For reconstructing a propagated action potential, one 
has to consider local circuit currents, which leads to [H: 

— (2) 

where R is the axon radius, p2 is the specific electrical 
resistance of the axoplasm, and x is distance along the 
axon. Combining eqs. (1) and (2), one obtains: 

R d^U ^ dU 



gNa{U~ENa)+gi{U~Ei). 



(3) 



Hodgkin and Huxley noted that this equation was not 
practicable to solve in this form. To simplify, they as- 
sumed that a steadily propagating solution exists, in 
which the shape of U against distance does not change 
in time, leading to: 

d'^U 1 d'^U 

(4) 



dx'^ 



where 9 is the velocity of conduction. From eqs. (3) and 
(4), one obtains: 
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dt 
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(5) 
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gNa{U-ENa)+9l{U-Ei). 

This is an ordinary differential equation and can be 
solved numerically if 9 is known. It is solved by guess- 
ing a value of 9, inserting it in eq. (5) and carrying out 
numerical integration starting from the resting state at 
the foot of the action potential. As the guessed 9 is too 
small or large, it is found that U approaches either +00 
or —00. A new value of 9 is then chosen and the proce- 
dure repeated until a correctly guessed 9 brings U back 
to the resting condition when the action potential is over 
Although the Hodgkin-Huxley model does lead to 
a propagating action potential, certain flaws were found 
in the formulation of the problem Q . Importantly, the 
model is not able to explain a number of mechanical, op- 
tical and thermodynamic phenomena accompanying the 
action potential, as well as it provides few clues for un- 
derstanding the Meyer-Overton rule 0, 



III. THE PRESSURE WAVE MODEL 

We suggest that the propagation of the nerve impulse is 
driven by a pressure pulse propagating in the axoplasm 
along the axon length. While several variations of the 
process can be envisioned, we suggest that it proceeds 
as shown in Fig. 2. A propagating axoplasmic pres- 
sure pulse distends the axon membrane, which causes 
mechanical activation of the membrane Na"*" ion chan- 
nels. Activated Na+ channels allow Na"*" ions to influx 
along the inward electrochemical gradient, thus locally 
depolarizing the membrane to the excitation voltage and 
initiating a local Hodgkin-Huxley voltage spike. The 
local transmembrane voltage spike activates membrane 
voltage-gated Ca^+ channels, leading to an influx of Ca^+ 
ions along their inward electrochemical gradient. The 
presence of free intracellular Ca^+ induces a contraction 
in actin-myosin filaments anchored to the membrane inte- 
rior. The filament contraction forces a radial contraction 
in cylindrical segments of the membrane (Fig. 3), which 
amplifies the propagating pressure pulse and thus com- 
pensates for its decay due to viscosity. (In section 5, it 
is shown that under physiological conditions axoplasmic 
pressure pulses decay over roughly 1 mm distances.) 

Figs. 4(a)-(b) illustrate a suggested process of me- 
chanical activation of Na"*" ion channels. The lipid bi- 
layer is stretched laterally as the propagating axoplas- 
mic pressure pulse causes a AR{t) increase in the mem- 
brane cylinder radius i?, corresponding to a AT(t) in- 
crease in the lipid bilayer tension r. Tension r can be 
expressed through the bilayer lateral pressure profile p{z) 

[131: T = f^^p{z)dz, where z is the depth within the 
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bilayer, and the integral is over the bilayer thickness 
h. The change AT(t) in the bilayer tension results in 
a change Ap(z, t) in the lateral pressure profile such that 

Ar(t) = Ap(z, t)dz. We suggest that these alter- 
2 

ations in the bilayer lateral pressures modulate confor- 



mational opening of Na+ channels imbedded in the bi- 
layer. Alternatively, it is plausible that Na+ channels are 
activated by direct mechanical links (such as filaments) 
tethering the ion channel structures to the interior cy- 
toskeleton that is perturbed by the propagating pressure 
pulse, as shown in Fig. 4(c). 



According to the model presented, an initial influx of 
Na+ ions through mechanically activated ion channels de- 
polarizes the membrane to the excitation voltage, leading 
to the development of a local Hodgkin-Huxley voltage 
spike. Further, the change in voltage across the mem- 
brane activates membrane voltage-gated Ca^+ channels, 
allowing entry of Ca'^^ ions into the cell, which initiates 
a contraction in the filament network anchored to the 
membrane interior. Ca^"*" is a good candidate for medi- 
ating the Hodgkin-Huxley voltage spike into the filament 
contraction for several reasons. First, Ca^+ ion channels 
are found in virtually every excitable cell [28j , and many 
types of Ca^+ ion channels are directly and rapidly gated 
by voltage [29| . Furthermore, an increase in free intracel- 
lular Ca^+ is often associated with initiation of motion in 
cells, from motility in freely moving cells and muscle con- 
traction to synaptic vesicle release at synapses fSsj. Also, 
the free intracellular Ca^"'' concentration is typically ex- 
tremely low, while it is much larger outside the cell. Ca^+ 
ions that influx during the action potential should in- 
stantly increase manyfold the free Ca^^ concentration in 
the vicinity of the filament network attached to the in- 
ner membrane surface; that, coupled with the ability of 
Ca^+ to quickly induce conformational changes in pro- 
teins (such as upon its binding to actin filaments in the 
actin-myosin muscle complex), could provide a fast con- 
tractile response necessary to amplify a propagating pres- 
sure pulse. In this scheme, the Hodgkin-Huxley sodium 
and potassium currents and the voltage spike simply pro- 
vide a means for activating voltage-gated Ca^"*" ion chan- 
nels. Alternatively, it is possible that amplification of 
the pressure pulse proceeds through an electromechani- 
cal coupling mechanism, such as voltage-induced mem- 
brane movement (soj resulting from the Hodgkin-Huxley 
voltage spike. Note that, since the radial bilayer contrac- 
tion follows the pressure pulse at the same velocity, even 
a small contraction yielding a small amplification will 
accumulate over distance and will create a larger ampli- 
fying effect. In this context, it should be noted that a 
pressure pulse propagating in an inviscid incompressible 
fiuid enclosed in a viscoelastic tube, dissipates energy and 
decays because the pressure exerted by the tube on the 
fiuid during tube radial expansion is greater than that 
during tube radial contraction; therefore, the total work 
done by the tube on the fiuid is negative. In axons, the 
situation may be the reverse - a larger pressure exerted 
by the bilayer during its contraction may result in an 
overall transmission of kinetic energy to the axoplasm. 
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Axoplasmic pressure pulse 
distends lipid bilayer 



Radial contraction of the lipid 
bilayer amplifies propagating 
axoplasmic pressure pulse 



Contracting actin-myosin filaments force 
a radial contraction in the lipid bilayer 



t 



Influxing Ca activates contraction 
in actin-myosin filaments anchored 
to the membrane interior 



Ca^* influxes along the inward 
electrochemical gradient 



Bilayer stretch opens 
mechanosensitive Na* 
channels 



Na* influxes along the 
electrochemical gradient 
and depolarizes membrane 
to the H-H excitation voltage 



Voltage-gated Na* channels are 
activated and a local voltage spike 
develops according to H-H mechanism 



The voltage spike across the membrane 
opens voltage-gated Ca^* channels 



FIG. 2: Diagram of the suggested process of propagation of the axoplasmic pressure pulse and the action potential. 



Radial bilayer 
contraction amplifies 
the pressure pulse 



Axoplasmic pressure 
pulse distends lipid bilayer 




Pressure 
pulse direction 



FIG. 3: Amplification of the axoplasmic pressure pulse by a 
radial contraction of the axon membrane (drawn not to scale). 



IV. MEYER-OVERTON RULE AND THE 
ACTION OF ANESTHETICS 

Within the framework of the nerve impulse as an ax- 
oplasmic pressure pulse, we naturally arrive at a quanti- 
tative explanation of the Meyer-Overton rule for inhaled 
anesthetics [3 113, [13] ■ Essentially, the rule states that 
anesthetic potency is determined by the bilayer concen- 
tration of anesthetic, independent of its molecular iden- 
tity [l^]. Therefore, the Meyer-Overton rule can be ex- 
plained if there is a property of the lipid bilayer that is 
essential for nerve impulse propagation and that, in addi- 
tion, depends only on the bilayer concentration of anes- 
thetic. As is shown below, one such property may be the 
bilayer area expansion modulus, which affects both the 
decay and the velocity of axoplasmic pressure pulses. It 
is also shown that anesthetics may inhibit the mechani- 
cal activation of membrane Na"*" channels by the travel- 
ing pressure pulse (and therefore suppress the conduction 
of the action potential) in a manner consistent with the 
Meyer-Overton rule. 



Suppose that an anesthetic molecule is dissolved within 
the lipid bilayer (anywhere from the headgroup to the 
hydrophobic region), and it weakens and destabilizes the 
bilayer structure nonspecifically, e.g., through creating 
a defect in the packing of lipid molecules and thereby 
nonspecifically disordering lipid chains (Fig. 5(a)). Such 
disordering interaction will generally decrease the bilayer 
elastic area expansion modulus in the vicinity of the 
lipid/anesthetic molecular complex (it is well established 
that bilayers with larger lipid chain disordering generally 
display lower elastic moduli and tensile strengths ,31]). 
Using a simple macroscopic model of the membrane j3l| , 
shown in Fig. 5(b), the total membrane area expansion 
modulus can be expressed as: 



K=(^ 
\Km 



K 



L/A 



(6) 



where the right-hand side of the equation is a combi- 
nation of the area expansion moduli of the two compo- 
nents, namely, non-anesthetized membrane Km (which 
describes the elasticity of the lipid bilayer, embedded 
proteins and other native structures) and a purported 
lipid/anesthetic molecular complex i^^i/yi {K[^/A < Km), 
scaled by their area fractions om and a^/A- Assuming 
that anesthetic molecules are dissolved only in the lipid 
segments of the membrane and that they do not cluster 
together, the area fractions om and a^/A are given by: 



flL/A ^ I - aM ^ TLA ■ h ■ A 



L/Ai 



(7) 



where ua represents the concentration of anesthetic 
molecules within the membrane, in units of molecules per 
volume, h the bilayer thickness, and Ai^/a the membrane 
area occupied by a single lipid/anesthetic molecular com- 
plex. Note that A^ /a and K^ /a are independent of the 
chemical identity of the anesthetic, due to the purported 
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Channel closed 





Change in membrane 
lateral pressures 
activates ion channei 



Channel open 




Channei open 

Pressure pulse 

FIG. 4: Mechanical activation of Na"^ 
traveling pressure pulse (see text). 



Channel is activated 
by fiiaments tethering 
channel structures 
to axoplasm interior 



ion channels by the 



nonspecificity of anesthetic/lipid interaction. Combining 
eqs. (6) and (7): 



K 



1 - UAhALjA^ ^ UAhA^iA 



K 



M 



K 



LI A 



(8) 



where the right-hand side of the equation does not de- 
pend on the identity of the anesthetic, but only on its bi- 
layer concentration n^- As will be shown in the next sec- 
tion, both the velocity and the decay length of axoplas- 
mic pressure pulses vary as \fK. Therefore, a decrease 
in the area expansion modulus K due to an increase in 
the anesthetic concentration ua will reduce the pressure 
pulse velocity and increase its decay, which, when ua is 
of sufficient magnitude, should lead to an inhibition of 
the pulse propagation. From the form of eq. (8), it fol- 
lows that anesthetics acting in such a way will comply 



Anesthetic molecule creates 
a defect in the lipid chain 
pacl<ing order 



Lipid bilayer 



Lipid bilayer 




Area expansion modulus is altered in the 
vicinity of an anesthetic molecule 



FIG. 5: (a) Disordering of lipid packing order by an anesthetic 
molecule dissolved in the bilayer. (b) Pockets of altered area 
expansion modulus in the vicinity of anesthetic molecules. 



with the Meyer-Overton rule. Prom the line of argument 

above, it is clear that any nonspecific anesthetic/lipid in- 
teraction that results in decreasing K with increasing ua 
and thus leads to a block in nerve impulse conduction 
will be consistent with the Meyer-Overton nde. 

Another possible mechanism of action of anesthetics 
consistent with the Meyer-Overton rule relies on inhibi- 
tion of mechanoactivation of Na"*" ion channels by the 
traveling pressure pulse, in cases where such activation 
is modulated by changes in the membrane lateral pres- 
sures (Fig. 4(b)). As before, we assume that anesthetic 
molecules dissolved in the bilayer weaken it nonspecif- 
ically and that increasing anesthetic concentration ua 
decreases the membrane area expansion modulus K. As 
will be discussed in the next section, the potential en- 
ergy of a propagating axoplasmic pressure pulse is mainly 
stored in the membrane strain. The membrane strain en- 
ergy, _Ep, is proportional to LAr^/if , where L is the pulse 
length and At is the peak increase in the membrane ten- 
sion r from equilibrium. Since L ~ \fK (assuming a 
fixed pulse duration, see eq. (14)), for the same pulse 
energy (Ep fixed), it follows that At ~ K^l^. Therefore, 
a decrease in the area expansion modulus K will result 
in a reduced tension increase At. At sufficiently high ua 
(low K), the decreased At may not suffice for activation 
of mechanosensitive Na+ channels, suppressing the pro- 
cess of pulse amplification and leading to a block of nerve 
impulse conduction. It should also be noted that, if the 
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anesthetic/lipid interaction changes the hpid bilayer area 
without affecting its area expansion modulus, the equi- 
librium tension in the membrane may change, which may 
affect gating of mechanosensitive Na"*" channels, possibly 
leading to conduction block. This would also conform 
to the Meyer-Overton rule if the anesthetic/lipid inter- 
action is nonspecific. In this context, it should be noted 
that an anesthetic/lipid interaction altering a mechani- 
cal property of the lipid bilayer in a nonspecific manner 
would probably involve a defect-like structure, in which 
physical properties of the anesthetic compound are not 
important provided that the compound satisfies certain 
criteria for creating the defect (e.g., being hydrophobic, 
nonpolar, and nonreactive and possessing certain spatial 
and mass dimensions). 



V. PROPAGATION OF AXOPLASMIC 
PRESSURE PULSES 

We consider the propagation of small amplitude, axi- 
ally symmetric pressure pulses in a viscous compressible 
axoplasm enclosed in a circular cylindrical thin-walled 
distensible membrane. The theoretical analysis of waves 
in viscous fluids enclosed in tubes has been presented 
in [H, m, [131 . For the limiting case of a viscous com- 
pressible axoplasm in a rigid tube, the pressure pulse is 
an axoplasmic density disturbance that propagates along 
the axon similarly to a sound wave packet, with the pulse 
potential energy stored in the axoplasm bulk deformation 
[32I [3^ . In the opposing limit of a viscous incompress- 
ible axoplasm in a distensible membrane, the pressure 
pulse manifests itself as an increase in the axon diameter 
and associated membrane area expansion, with t he p ulse 
potential energy stored in the membrane strain [sl, • 
Let us consider a pressure pulse of central frequency lo 
propagating in the axoplasm of density p, compressibil- 
ity K and dynamic viscosity /z, enclosed in a cylindrical 
membrane of radius R, thickness h and area expansion 
modulus K. Let us also assume that for the deformation 
of area expansion (when the membrane material behaves 
in an elastic solid manner), the membrane Young's mod- 
ulus is E and the Poisson's ratio is v. If viscosity is 
neglected, the pulse propagates with the velocity [35| : 



1 



Wo = 



(9) 



For typical unmyelinated axons, i.e., assuming K — 
0.8 N/m (twice the single plasma membrane area ex- 
pansion modulus of 0.4 N/m to account for the 
combination of the axon and adjoining glial cell mem- 
branes), K = 4.04 • 10-1° Pa-i (water at 38 °C [13), 
p — 1050 kg/m'^ (axoplasm [11]), and 2R — 1 /im, eq. 
(9) yields the pressure pulse velocity of 28 m/s. Note 
that if the axon membrane is assumed to be indistensible 
(^ ^ k), eq. (9) yields a pulse velocity of 1535 m/s. 



A. Viscous axoplasm 

Viscous forces in the axoplasm and in the axon mem- 
brane in general decrease the velocity of a xop lasmic pres- 
sure waves and introduce a decay [12, [H, [lj| ■ The mem- 
brane viscous tension is characterized by the ratio of the 
surface viscosity (on the order of 10~^ N-s/m [HI) and 
the time scale of the pressure pulse (here assumed to be 
close to the action potential duration, about 1 ms), yield- 
ing 10"'^ N/m. Since this value is much smaller than the 
membrane elastic area expansion modulus, 0.4 N/m, the 
effects of membrane viscosity on the pulse propagation 
can be neglected [1^. Concerning the viscosity of the 
axoplasm, for the general case of a viscous compressible 
fluid enclosed in an elastic tube, the dispersion equation 
is quite complicated and has not been solved in closed 
form. However, the equation is simplified and can be 
solved in the high viscosity limit, defined by: 



(10) 



Here we evaluate eq. (10) for typical axons. As before, 
we assume that the duration of the axoplasmic pressure 
pulse is similar to that of the action potential, about 
1 ms, yielding to « 3142 rad/s. Axoplasm viscosity /i, 
as is relevant for perturbations associated with small- 
amplitude axoplasmic pressure pulses, is assumed to be 
similar to the viscosity of water, /x = 6.82 • 10~^ Pa-s at 38 
°C [13 • Taking 2i? = 10/im (the typical internal myelin 
diameter for Aa myelinated fibers) and p = 1050 kg/m'^ 
yields a = 0.35, which should be small enough for high- 
viscosity approximation to the dispersion equation to 
hold [34] for axons of similar or smaller diameters. When 
a ^ 1 (2i? ^ 29/im), the effects of viscosity are small, 
and the pressure pulse propagates with little decay and 
with the velocity described by eq. (9). 

When eq. (10) holds, the phase velocity of axoplasmic 
pressure waves (the velocity of the harmonic waves) is 
given by i32„,33„^]: 



(11) 



where a and vq are given by eqs. (10) and (9), respec- 
tively; c — 1 when ^ ^ k (the limit of a rigid mem- 
brane), while c — when ^ ^ k (the limit of a 
soft membrane). The factor c accounts for movement of 
the membrane in the axial direction caused by the viscous 
axoplasm in the limit of a soft membrane; here and below 
we assume = 0.5 for the incompressible lipid bilayer, 
leading to ~ 1.15. The velocity of the pressure 

pulse (the "group" velocity) is given by Vgr = where 
k is the wavenumber. From eqs. (9) - (11), and using the 
relations Vph = j and Vgr = we obtain the velocity of 
the axoplasmic pressure pulse in the high viscosity limit 
(10): 



Ugr 



— c ■ a ■ vo = c ■ R 



(12) 
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Under the same limit, the distance over which the pulse 
amplitude decreases e-fold (the decay length) , is [12, [H, 




(13) 



B. Myelinated axons 



In the Hodgkin-Huxley theory of propagation of the ac- 
tion potential, the role of the myelin sheath is in insulat- 
ing the membrane electrically and in decreasing the mem- 
brane electrical capacitance; both contribute to more 
rapid conduction of the action potential. However, the 
myelin sheath is also quite rigid. About 80% of dry 
myelin content is lipid, with a roughly 2:2:1 molar ra- 
tio of the three major lipids, cholesterol, phospholipids 
and galactolipids [4^. Cholesterol is known to increase 
the strength and elastic modulus of lipid bilayers, with 
the highest modulus and highest strength achieved for 
50 mol% cholesterol 



31|, i.e., rather close to the con- 



centration in myelin. The remaining 20% of dry myelin 
content consists of protein, including various types of 
fibrous proteins that can form rigid structures. Alto- 
gether, for the deformation of circumferential strain, the 
myelin elastic area expansion modulus is about 2 N/m 
for a single myelin layer of 4 nm thickness [3l[. Noting 
that K = Eh (for v ~ 0.5), the myelin effective Young's 
modulus is i? = 5 • 10^ Pa. The ratio of the axon di- 
ameter (internal myelin diameter) to the nerve fiber di- 
ameter (external myelin diameter) is typically about 0.7 
[4]| . implying a myelin thickness of « 0.43i?. There- 
fore, ^ = -H « 9.3 • 10-9 Pa-i. This is about 23 
times greater than the compressibility of water at 38° C, 
K — 4.04 • 10^^" Pa"^. Hence, vq, as used in eq. (12), is 
only -\/24 « 4.9 times less than the speed of sound in wa- 
ter (see eq. (9)). In other words, the myelin sheath is so 
rigid that its distensibility decreases the pulse speed only 
about 4.9 times compared to the theoretical maximum if 
axon walls were absolutely rigid. For typical myelinated 
axons, e.g., cat myelinated axons of 10 /im fiber diameter 
at 38 °C, eq. (12) yields Vgr = 88 m/s (using S = 5 • 10® 
Pa, 2R = 7fim, h = 0A3R, v = 0.5, w = 3142 rad/s, 
IJL = 6.82 • 10^4 Pa-s and k = 4.04 • 10"^° Pa^^). Given 
the roughness of the estimate, this theoretical value of 88 
m/s is in excellent agreement with the measured action 
potential velocity of 60 m/s for these fibers f45| . From eq. 
(13), the decay length for the propagation of the pulses 
in myelinated segments of such a fiber (where the pulse 
propagates passively) is 14 mm. Assuming that the dis- 
tance between the nodes of Ranvier is 1 mm (100 times 
the external diameter of myelin), the pulse amplitude de- 
cays 7% over the internodal distance. 



C. Unmyelinated axons 

Unmyelinated nerve fibers usually have diameters of 
0.1 - 1.2 /im. The action potential conduction velocity, 
in m/s, is given approximately by v « 1800-\/7?, where 
R is the axon radius in meters [43| (values quoted in the 
literature range from V « lOOOv^ H to w « 3000/R 
[4^). Therefore, for an axon of 1 p,m diameter, conduc- 
tion velocity is about 1.27 m/s. Assuming 2R = 1 /im, 
K = 0.8 N/m (twice the single plasma membrane value 
of 0.4 N/m l3Q]), and lo, 11, n and v as before for myeli- 
nated axons, eq. (12) yields a pressure pulse velocity of 
Vgr = 1.11 m/s, very close to the experimental value. 
From eq. (13), the decay length for this fiber is 0.18 mm. 
It is truly remarkable that, for both myelinated and un- 
myelinated fibers, the predicted velocity of propagation 
of axoplasmic pressure pulses is very close to the mea- 
sured velocities of nerve impulses. 



D. Dependence of propagation velocity on fiber 
diameter 



As can be seen from the above numerical estimates, for 
both myelinated and unmyelinated axons the condition 
^ ^ K is true; and therefore, eq. (12) simplifies to: 



l2Ru;K 
3/i 



(14) 



This relationship shows that for unmyelinated axons, the 
axoplasmic pressure pulse velocity is proportional to the 
half-power of the axon diameter, assuming that the pulse 
duration, the membrane area expansion modulus and the 
axoplasm viscosity do not vary with axon diameter. For 
myelinated axons, introducing 7 = /i/i? (the ratio of the 
myelin sheath thickness to the axon radius) and using 
K = Eh, we can rewrite eq. (14) as: 



R\ 



l2ujEj 
3fi 



(15) 



Therefore, for myelinated axons, assuming that pulse du- 
ration and specific properties of axoplasm and myelin 
are independent of diameter and that the myelin sheath 
thickness scales in proportion to the fiber diameter (i.e., 
7 is constant), the pressure pulse velocity scales linearly 
with the fiber diameter. It is interesting that a linear 
dependence of pressure pulse velocity on fiber diameter 
is also obtained for myelinated axons without assuming 
K < from eq. (12): 



R. 



(16) 



In summary, given that specific properties of axoplasm 
and myelin are independent of fiber diameter and that 
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the myehn sheath thickness scales in proportion to fiber 
diameter, the axoplasmic pressure pulse velocity scales 
as Vi? for unmyelinated axons and as R for myelinated 
axons, similar to what is measured experimentally [43!, 

m. 

E. Dependence of propagation velocity on 
temperature 

The duration of the action potential increases with de- 
creasing temperature, e.g., with Qio (the ratio of the 
velocity at one temperature to the value at a tempera- 
ture 10 °C lower) of 3.4 at 37 °C for cat vagus myelinated 
fibers [i^ . Assuming that the duration of the axoplasmic 
pressure pulse is similar to that of the action potential, 
taking the change in /i from 37 °C to 27 °C as for water, 
Qio = 0.81 [131, and assuming that the other parameters 
in eq. (12) do not change with temperature, we obtain 
the predicted Qio of 2 for the nerve impulse propagation 
velocity. This is rather close to the experimental value of 

1.6 [ig. 

VI. DISCUSSION 

It has recently been suggested [1, [1] that the Hodgkin- 
Huxley model [1| of the action potential does not provide 
a satisfactory description of the nerve impulse because it 
does not include the mechanical d, 0, M, H, and op- 
tical [H, [l3| changes associated with the action poten- 
tial, as well as it is inconsistent with observed reversible 
changes in temperature and heat [ll|, [l^, fisl. [13 . fisl . Fur- 
thermore, the famous Meyer-Overton rule [l8l. Il9l [20| is 
seemingly inconsistent with the Hodgkin-Huxley theory, 
as it implies an action of anesthetics other than binding 
to ion channels The authors suggested a theory of 
the action potential based on solitons propagating in the 
axon membrane [1, |^ . Although this theory is in princi- 
ple able to account for the abovementioned phenomena, 
apparently it cannot explain changes in the propagation 
velocity that are related to axon diameter. 



Following earlier work [2, l3|, we have here proposed 
an alternative theory of the propagation of the nerve im- 
pulse. We suggested that the traveling transmembrane 
voltage spike (the action potential) is preceded by a pres- 
sure pulse propagating in the axoplasm. We have shown 
that, under physiological conditions, such pulses decay 
over roughly 1 mm distances; therefore, sustained prop- 
agation requires a mechanism of amplification. Based on 
known axon properties, we have proposed several am- 
plification mechanisms. Each of these mechanisms in- 
volves the Hodgkin-Huxley voltage spike as a mediator, 
thereby providing the spike a functional role within the 
model. Due to the presence of the axoplasmic pressure 
pulse, the model is consistent with the mechanical, opti- 
cal and thermal ch ang es associated with the action po- 
tential [i, 0, i i, [M [nl, [il [il, Q [il, [M [i3- im- 
portantly, the simple and powerful Meyer-Overton rule 
[Tsl [TqI . |20| . valid over 5 orders of magnitude, currently 
appears to be difficult, if not impossible, to interpret 
within the Hodgkin-Huxley theory. Within the presented 
framework, however, the Meyer-Overton rule is readily 
explained based on properties of lipid bilayers. A simple 
test of the presented theory is its ability to predict the 
velocity of the nerve impulse. We have shown that, over 
2 orders of magnitude (for myelinated and unmyelinated 
axons), the predicted velocity is close to that which is 
experimentally observed, despite the fact that the preci- 
sion of calculations is limited by current knowledge of the 
mechanical properties of axons and of the duration of the 
axoplasmic pressure pulse. The dependence of nerve im- 
pulse velocity on fiber diameter and temperature is also 
reproduced well by the model. A transmembrane voltage 
spike such as the action potential, however, is not a neces- 
sary feature of the model; therefore, it is suggested that 
some neurons may not exhibit the voltage spike, while 
still transmitting the nerve signal over short distances. 
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